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Abstract 

We generalize the Lee-Suzuki iteration method for summing the folded di- 
agram series to the case where the unperturbed model-space energies are 
non-degenerate. A condition is derived for the convergence of the iteration 
scheme and this depends on the choice of the model space projection oper- 
ators. Two choices are examined, in the first the projection operators are 
defined in terms of the unperturbed states and in the second they are defined 
in terms of the eigenfunctions obtained at each stage of the iteration. As is 
illustrated by calculations with a simple model, the second procedure gives 
the better convergence and, by suitable choice of the starting energies, allows 
the reproduction of any subset of the exact eigenvalues. 

PACS number: 21.60.Cs 



1 Introduction 



In nuclear, atomic and chemical physics it is usually necessary to recast the 
full many-body problem in the form of an effective interaction acting within a 
chosen model space for which the eigenvalues can be obtained exactly. Much 
work has been carried out on this topic, both as regards formal questions 
and actual calculations, see refs. [1-5] for the nuclear case. The formalism 
consists of a completely linked perturbation series which contains both non- 
folded and folded diagrams. For a given set of non-folded diagrams, the 
folded diagram series can be summed by using either the Krenciglowa-Kuo 
(KK) technique || or the Lee-Suzuki (LS) method |7|]. Both these methods 
employ completely degenerate unperturbed energies for the model space, 
however their convergence properties are different. The KK approach, when 
convergent, yields the eigenvalues for those states which have the largest 
overlap with the chosen model space. On the other hand the LS method 
reproduces those eigenvalues which lie closest to the chosen unperturbed 
energy. 

Since actual single particle energies in, for example the (sd) or the (pf) 
shells, are far from degenerate, it is clearly desirable to use a formalism which 
is not restricted to exact degeneracy. This would allow one to treat the one- 
body terms as unperturbed energies rather than introducing artificial energy 
shifts so as to rewrite the problem in degenerate form. Further the KK and LS 
methods yield only certain of the exact eigenvalues which is not, in general, 
desirable. This restriction is not present in the exact representation of the 
complete many-body problem as a series of non-folded and folded diagrams 
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since this contains information regarding all of the true eigenvalues. Given a 
model space of dimension d it should be possible to obtain any selection of 
d eigenvalues from the complete set of true eigenvalues. We shall show that 
this is indeed possible using the formalism developed in this paper which is 
expressly designed for non-degenerate unperturbed energies. 

The organization of this paper is as follows. In sec. 2 we briefly outline 
the standard LS formalism and establish notation. The generalization of the 
LS approach to the non-degenerate case is discussed in sec. 3. Since the 
solution of the equations requires iteration, we give in sec. 4 a criterion for 
the convergence of the iteration. In order to assess the present approach and 
to compare with the KK and LS methods, we need to study a case where 
the exact results are known, i.e. a model. This is the subject of sec. 5. Our 
concluding remarks are given in sec. 6. 

2 Outline of the Lee-Suzuki Method 

In the usual way we write the full Hamiltonian, H = Hq + V, where Hq is the 
unperturbed Hamiltonian and V is the perturbation. We use a basis in which 
H Q is diagonal and define an operator P which projects the Hilbert space 
onto the chosen model space of dimension d. The complementary operator 
Q projects onto the remainder of the Hilbert space, thus P + Q = 1, PQ = 0. 
Then it is straightforward to obtain formal expressions for the standard non- 
Hermitian effective interaction, which we denote here by R. It can be derived 
by a number of different methods ^ ^ § and can be written in various forms. 
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Here it is convenient to use 

R = PVP + PVQu , (1) 
where the operator oj obeys the equation 

QVP + QHQu - uPHP - uPVQu = . (2) 

Then the model-space eigenvalue equation, yielding d of the true eigenvalues 
labelled E p , can be written 

{PH G P + R)\4> P )=E P \4> P ) , (3) 

where the model space wave function is the projection of the true wave func- 
tion on the model space, i.e., <p p = P^ p . 

In this section we consider a system with degenerate unperturbed ener- 
gies, thus 

PH P = e P . (4) 
In this case the equation for u becomes 

(e - QHQ)u + uj(PVP + PVQuj) = QVP , (5) 
which can be rewritten in the form 

Substituting eq. (6) into eq. (1), we have 



R 



1 

l + PVQ — — co 

e - QH Q 



Q(eo) , (7) 
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where we have defined the Q-box according to 

1 



Q(e ) = PVP + PVQ- 



-QVP . 



e - QHQ ~* ' * ' (8) 
The LS method J7j generates the solution of eqs. (6) and (7) by iteration. 
The n th iteration is obtained from the (n — l) st iteration by the following 
equations 

1 



Rn : 

and 



P + PVQ 



e - QHQ 



1 



1 



Q(eo) , 



■ ^n-\Rn ■ 



e -QHQ^ eo-QHQ 
Defining Uq = 0, the n th iterative solution, (n > 2), is given by 



(9) 



(10) 



Rn= P-Ql 



QlRn-X — ■ ■ ■ - 
ra—l n— 1 



^ " Ql ~ E <9m II R k 

m=2 k=n—m+l 



Qn-lR2R3 
1 

Q , 



Rn-l Q 



(11) 



where Q m for m = 1, 2, ... is given by the m derivative of the Q-box, namely 

«- = <- irpv '° (^ef 1 = ' (12) 

The sequence {R\, R2, . . .} is generally convergent (see sec. 4) so the solution 
for the effective interaction corresponds to R = Roq. The rate of convergence 
and, indeed, which of the exact eigenvalues are obtained from R^ will de- 
pend on the unperturbed energy eo- It is important to realize that this is a 
parameter at our disposal. Thus the Hamiltonian and eq. (5) are unchanged 
if we shift eo to eo + e' and compensate for this with a corresponding shift of 
the perturbation V to V — e'P; the quantity e' is clearly arbitrary and can 
be used to optimize the convergence of the iteration procedure. 



3 The Non-degenerate Case 



3.1 General Iterative Solution for Effective 
Interactions 

Let us define projection operators, P a , which act in the model space and are 
such that 

P = £P Q and P a P p = 5 al3 P a . (13) 

a 

It follows that 

P a P = PP a = P a . (14) 

We can then discuss the general situation where the model-space eigenvalues 
of Hq are not completely degenerate. As in the degenerate case, we have 
the freedom to modify the unperturbed P-space Hamiltonian arbitrarily and 
make a compensating change in the P-space part of the perturbation. Such 
shifts will affect the convergence of the iteration and will determine which of 
the exact eigenvalues are finally obtained. Thus we choose 

ph' p = £* p * and py'p = py p + PHoP - e « p « • ( 15 ) 

a a 

This does not exclude the possibility that some degeneracy is still present, 
indeed the formalism may be applied to the completely degenerate case in 
which case one regains the standard LS method of the preceding section. 

A general perturbative expansion has been given in ref. for the non- 
degenerate case within the framework of the Q-box formalism. Here we wish 
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to generalize the iterative scheme of sec. 2. To that end we substitute PH' Q P 
in eq. (15) into eq. (2) and multiply by P a from the right, yielding 

QVP a + QHQuP a - e a uP a - u(PV'P + PVQu)P a = , (16) 



from which it follows that 
1 



uiP a 



QVP a - 



■ UjRP n 



(17) 



e a - QHQ " e a - QHQ 
Using eq. (13) and noting that ujP = ui we have a generalization of eq. (6) 
to the non-degenerate case, namely 



UJ 



uRP a 



(18) 



^e a -QHQ^ „ £«,- QHQ 

Multiplying eq. (1) on the right by P a , using eq. (17) and noting that 
R = RP = J2 a PPa, we have a formal solution for R given by 

1 



1 + PVQ- 



■ UJ 



e a - QHQ 
where Q(e a ) is similar to eq. (8), namely 



Qi^a)Pa i 



Q(e a ) = PV'P + PVQ- 



1 



:QVP 



(19) 



(20) 



t a - QHQ 

We now set up iterative equations for uj and R. In the most general case 
the projection operators may vary according to the iteration, i.e. P a — > P™. 
This means that Q — > Q n since it is dependant on the iteration through the 
quantity PV'P = PVP + PH P - £ a e Q P™. We can write the iterative 
equations in the form 

1 



Rn~Yl 

a 

and 



i -i 



P + PVQ 



e a - QHQ 



Q 



n—l, 



e a)P a ) 



QHQ 



QVP™- 1 - 



QHQ 



^n-lRnPa 



n—l 



(21) 



(22) 



This is just the generalization of the LS eqs. (9) and (10) to the non- 
degenerate case. We define P° in terms of the basis states of the original 
unperturbed Hamiltonian, thus 



P° = |a;) (a | where H 
Taking u = 0, we have the sequence 

i2i = E0°(O^, 



a) = e n \a) 



R2 = E 



P 



-P _ e /3)"P/3 ~~ E<M e a, e l)Pp^P^ 

P Pi 



a 

and in general i? n , for n > 2, is given by 



(23) 

(24) 
(25) 

-i 

Q 2 (e a )Pl(26) 



P — E Ql( € a' 6 p)Pp 2 ~ E Qz( 6 on e /3, e ^i)Pp 3 Rn-lP- 



n-2 

1 



P 



Pi 



-i -1 



— '^2Qn-l(ta, £/3, e 7) ' ' ' e A, e M )Pgi?2-P^-R3 ' ' " -Rn-2-P" ^ Rn-lP^ 

n-1 



n-2 



= E 



^-EQi( e «' e /?) ? (3 2 ~E E ^( f «.eA,--,e/)J 



pn-m-1 

X 



Here we have defined 



Q 1 (e 1 ,e 2 ) = -PVQ 



m=2 fcfo-Pm 

k=n— m+1 



QVP 



(27) 



Cl - QiJQ e 2 - Q#Q 
fifeiizflisi if ei ^ e 2 
if ei = e 2 



ei— £2 



(28) 
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and, in general, 

$-(«.,«., • ■ •,«-.) = (-D-PVO— 1^ ■ ■ • 

e m +i _ QHQ^ ^ ^ 

3.2 Choice of Projection Operator 

The simplest choice is to keep P a fixed at its initial value, namely 

PS = P°=\a){a\. (30) 

This means that PV'P = PVP so that Q n is independant of n and is the 
same as Q of eq. (8). We refer to this as the generalized Lee-Suzuki approach 
(GLS). We can easily regain the results of sec. 2 in the completely degenerate 
case where e a = eo for all a. Thus Q m (ei, 62, • • • , e m+ i) becomes simply the 
Q m of eq. (12). This means that the summations in eq. (27) are just of the 
form Y^3i P/3i — P anc i since Q m P = Q m w e return to eq. (11). 

The other choice we consider for the projection operators is to write 
them in terms of the model space wave functions obtained at each iteration. 
Specifically our n th approximation to the exact eq. (3) is 

(E^r 1 + ^|0^)=^|0^), (31) 

where E% is the n th approximation to the true eigenvalue E a . The vectors 
being simply the projections of the true wave functions onto the model 
space, are not orthogonal. However it is well known that their biorthogonal 
complements can be defined such that = 5 a p- We then define 

CHCXCI, (32) 



and it is easily verified that eqs. (13) and (14) are satisfied. It is necessary 
to specify which e a is associated with a given P a in eq. (31). This we do in 
the obvious way by ordering the unperturbed energies e x < e 2 < • • • < and 
the approximate energies E™ < E% < ■ ■ ■ < and making the association 
with the corresponding eigenvectors. We refer to this as the generalized 
Lee-Suzuki method with self-consistent basis (SCGLS). In the completely 
degenerate case J2 a e aP% becomes e J2 a P£ = e oP an d the approach reverts 
to the standard LS method by the arguments given before. 

3.3 Evaluation of Q m 

The operator Q m (ei, e 2 , • • • , e m +i) is the basic element needed to construct 
the effective interaction for the non-degenerate case and it can be expressed 
as a linear combination of the standard Q-boxes according to 

m+1 

Qm(ci, £2, • • ■ , Wi) — X! Ck(ei, e 2 , • • • , e m +i)Q(ck) , (33) 

fc=i 

where 

Cfe(ei, £2, • • • , e m +i) = (efe — ei)(efe — e 2 ) • ■ • (e^ — €k-i){tk — Cfe+i) • • • 

• • • (e k - e m+ i) . (34) 

For the derivation of eq. (33) we have used the equality 

\ei-QHQ) ■ ■ ■ (e^-QHQ)]" 1 = ^(-i^ ^h^l^A ^) 

which is easily proved by induction. It should be noted that the term PV'P 
in the Q-hox. (eq. (20)) gives no contribution to eq. (33) because of the 
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relation 

m+l 

c k(ei,e 2 , • ■ ■ , e m +i) = , (36) 

fc=i 

which can be proved directly from eq. (34). 

So far we have implicitly assumed that the energies (ei, e 2 , ■ ■ ■ , e m +i) are 
all different. However in carrying out the summations we shall need to con- 
sider the case where two or more of the energies refer to the same state 
and there may also be some degeneracy present. Suppose, for instance, that 
ei = e 2 so that the formal limit 6\ — > e 2 leads to the derivative as in eq. 
(28). This is actually evaluated by calculating Oi( e i> e i + 6) or > m general, 
Q m (ci, £i + 6, e 3, • • • , Wi)) where 5 is small in comparison to e\. Similarly if 
ei = e 2 = £3, we would use Q m (ei — 5, e±, €1 + 6,64, ••• , e m +i), and so on. Thus 
the fundamental quantity we need is Q(ek +p6), where p — 0, ±1, ±2, . . . and 
this only differs from the degenerate case in that here more than one value of 
€k needs to be considered. Thus in applying eqs. (33) and (34) any degener- 
acy should be broken by a small amount so that the appropriate derivative 
is implicity evaluated. 



4 Convergence of the Iterative Solution 

We now discuss the convergence condition for the generalized LS iteration 
scheme. Let 5u n ^i and 5uj n be the deviations from the exact solution u of 
eq. (18) in the (n — l) st and n th iterations respectively. Thus, by definition, 

u n ^i=u + 6u n _i , (37) 
uo n = uj + 5uo n . (38) 
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We first discuss the case where the projection operator is fixed, i.e., P™ = P°. 
Now setting R n = R + 5R n , we find from eq. (1) that 5R n = PVQ5u n . Then 
to first order in 6u, using eq. (22), we have 

[e Q - (QHQ - luPVQ)] 5u n P° = -5u n ^RP* . (39) 

Now the term (QHQ — ujPVQ) is the Q-space effective Hamiltonian 0. The 
eigenvalue equation for this case can be written in the form 

(ip q \(QHQ - uPVQ) = E q (^ q \ . (40) 

The eigenvalues E q above and E p in eq. (3) agree with two of the exact 
eigenvalues of the Hamiltonian H. Now dividing eq. (39) by the operator in 
square brackets on the left and multiplying by (ip q \ on the left and \<f> p ) on 
the right, we obtain 

(^l<w n p°l0 P > = (vg<W-iPP> P >(p g - e a y} (4i) 

This can be simplified by using \(f)), the biorthogonal complement to \<p), and 
inserting J2 P ' |0p'}(<?Vl = 1 since 

(4> p ,\RP° = (E p ,-e a )(4> p ,\P°. (42) 

Performing these manipulations and summing over a in eq. (|4T|), we finally 
obtain 

(ip q \5u n \<f> p ) = ^(^|5^ n _i|0p/)Z*, p , (43) 
P > 

where we have defined 
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It is useful to write eq. ( f4~3| ) in an obvious matrix notation as 

5uj n = 5w n _i ■ Z . (45) 

Then in order that the iteration be convergent it is necessary that the norm 
|| 5(jJ n || be smaller than || 5u? n _i || for n greater than some integer N. This 
means that a sufficient condition for convergence is 

||Z||<1. (46) 

Here we employ the Hilbert, or spectrum, norm which is defined for an ar- 
bitary matrix X as || X ||= maxv^Ai, where Aj are the eigenvalues of the 
matrix X^X. In the completely degenerate case we see that eqs. (f44|) and 
( f4~6"|) reduce immediately to the LS convergence condition 
e — E p 



<1. (47) 

Thus eq. ( f4~6| ) represents the generalization of the LS condition to the case 
of a non-degenerate P-space. 

The convergence that is obtained will depend on the choice of projection 
operator. In the GLS case the condition (46) may be quite complicated, since 
eq. (44) involves (</y |-P°|0p) = |a)(a| </> p ). However, if the states that 



we wish to obtain, \<p p ), lie largely within the model space the non-Hermiticity 
of the effective interaction will be small (see the discussion in ref. |1(|) and 



\(p p ) will not differ greatly from \(j) p ). Further if a given state p contains a 
large component of \a) we can expect Z to be close to a diagonal matrix. 
Then the convergence condition (|46| ) should be satisfied if the unperturbed 
energy e a is close to E p , but distant from the Q-space eigenvalues E q . 
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We next discuss the convergence condition for the SCGLS approach. In 
this case we have to consider the deviation of the projection operator 
defined as 

5 pn = pn _ poo ^ (4g) 

with 

P? = \4>a)(ja\, (49) 

where \<f> a ) is the exact eigenstate in eq. (3). Since R n in eq. (21) can be 
written as 

R n = PV'^-VP + PVQoo n 

= PVP + PH P - ]T e a P2~ l + PVQu n , (50) 

a 

the deviation SR n is given by 

5R n = PVQ5uj n - ]T tJP^ 1 . (51) 

a 

Substituting eqs. (37), (38) and (51) into eq. (22), we see that, to first order, 
the terms with the deviation are all canceled and we have 



a 



[e a - {QHQ - uPVQ)]5u n P™ = -Su^RP™ , (52) 

which is just the equation obtained by replacing P° by P^ in eq. (39). 
Directly from eq. (52) or by replacing P° by P^ in eq. (44), we obtain the 
matrix Z^ p , as 

z * = fcf ) ' (53) 
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We finally may say that in order that || Z ||< 1 the following condition must 
be satisfied: 



max 



E p 



£ P — Eg 



<1, p=l,2,...,d, q = l,2,... . (54) 



The above condition is satisfied if and only if the d eigenvalues E p correspond 
to those true eigenvalues which lie nearest to the unperturbed energies e p with 
p= l,2,...,d. 

The magnitude of | |Z| | is equal to unity at points where e p = \{E P + E q ), 
i.e., one of the unperturbed energies is exactly halfway between a P-space 
and a Q-space eigenvalue. Apart from these points eq. (54) is obeyed, so that 
this iteration scheme will converge to those eigenstates whose eigenvalues lie 
nearest to the unperturbed energies e p . 



5 Test Calculations 

In order to obtain some assessment of the GLS and SCGLS methods, as well 
as to compare with the KK and LS methods, we need to study a model for 
which exact results are readily obtained. The model we shall use is a slightly 
modified version of one which was introduced many years ago by Hoffmann 



et al. [lTJ in order to study the intruder state problem. The Hamiltonian 



is taken to be H = H Q + V, where the unperturbed Hamiltonian matrix is 
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Hq = diag(l, 1,3,9) and the perturbation V is given by 



V = 



( 





5x 


—5x 


5x \ 




5x 


25x 


5x 


-8x 




—5x 


5x 


—5x 


X 




5x 


-8x 


X 


—5x J 



where x is a parameter that we shall vary. Obviously the eigenvalue problem 
for this model Hamiltonian can easily be solved. 

We shall take the lowest two states of H to be our model space. Our 
main concern is then to see whether an effective interaction acting in this 
model space is able to reproduce any pair of exact eigenvalues. In tables 1 
and 2 we show results obtained for x = 0.05 and 0.20. Physically these cases 
differ in that for x = 0.05 the model space states dominate the eigenvectors of 
the first and second states, whereas for x = 0.20 they dominate the first and 
third eigenvectors, i.e., increasing the value of x causes a "level crossing" to 
take place. Using the SCGLS method we show results obtained with various 
choices for the e a of eq. (15). We see that as the unperturbed energies change 
the solutions converge to different pairs of eigenvalues. As predicted in sec. 
4, the SCGLS yields the eigenvalues which are nearest to the unperturbed 
energies. The convergence rates are reasonable in all cases, with accuracy to 
four decimal places obtained in at most 12 iterations and usually many fewer 
are required. By adjusting the e a we are able to reproduce any pair of exact 
eigenvalues, regardless of the order of the eigenvalues or the magnitude of 
the P-space overlap of the eigenvectors. 

We have also carried out calculations with the GLS method where the 
projection operators are defined in terms of the unperturbed eigenf unctions 
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and do not vary with iteration. For x = 0.20 we were able to reproduce 
all pairs of eigenvalues as with the SCGLS, however for x = 0.05 we were 
not able to obtain convergence for three of the combinations. These were 
(Ei, P 4 ), (E 2 , E 4 ) and (E 3 , P 4 ). All of them involve state 4 which has a very 
small overlap with the model space and the quantity Z in eq. (44) depends on 
the overlaps ((f) p >\a) and (a\(fi p ) as well as the energy differences (e a — E p >) and 
(e a — E q ). Therefore it is possible to have || Z ||> 1, implying no convergence, 
even though |(e a — E p >)/(e a — E q )\ < 1. 

We now turn to a comparison of the different iteration methods- KK, LS, 
GLS, SCGLS. Briefly in the KK method || the effective interaction in the 
n th step is given by 

fin=EQKHC 1 )ri, (55) 

k 

where 

(ph p + Rn-.w,- 1 ) = ^rvr 1 ) • (se) 

This iterative process sums the folded diagrams to all orders and, if it is 
convergent, the states with maximium P-space overlap are obtained. For x = 
0.05 the two lowest states are the ones with maximum P-space overlap and for 
this case we can compare the four methods. In order to make the comparison 
we introduce a measure of the deviation of the calculated eigenvalues from 
the exact results, 

1 

(57) 

We plot A n versus the number of iterations n in fig. 1. Here we have taken 
€\ = 1.0 and e 2 = 2.0 for the GLS and SCGLS and a degenerate energy of 
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L k 



TP exact \ A 



eo = 1.5 for LS and e = 1.0 for KK. Of course all methods converge, but the 
SCGLS and KK techniques show the fastest convergence. Thus it appears 
that the use of a self-consistent projection operator yields the most rapid 
convergence. As another example, we take x = 0.20 and reproduce the first 
and third eigenvalues, which are the ones with maximum P-space overlap, so 
that the KK, GLS and SCGLS methods can be compared. Taking eo = 1.0 
for KK and e\ = 0.0, e 2 = 5.5 for GLS and SCGLS, we obtain the results 
in fig. 2. The convergence of the SCGLS method is markedly better than 
for the other cases. In particular the GLS is rather slowly convergent and 
shows an oscillatory behavior, although it does ultimately yield an accurate 
answer. In fig. 3 we make a comparison between LS, GLS and SCGLS at this 
value of x (using eo = 4.0 for LS and e\ = 2.5, e 2 = 5.5 for the other cases). 
These methods converge to the second and third eigenvalues because these 
lie closest to the unperturbed energies. Here the GLS and SCGLS show a 
better rate of convergence than the standard LS method, with the SCGLS 
giving the most rapid rate. 

Finally it is of interest to examine the starting energy dependance of the 
SCGLS method. We show in fig. 4 the effect of varying e 2 while keeping 
€i = 0.0 for x = 0.20. The calculations then yield the ground-state energy 
and either the second or third or fourth eigenvalue depending on e 2 . The 
solid curve in fig. 4 gives the minimum number of iterations, n min , required 
for A n < 10~ 4 . The convergence is good except at the points where e 2 = 
l(E 2 + E 3 ) = 4.11 and e 2 = |(P 3 + P 4 ) = 7.29, since here || Z ||= 1. This is in 
agreement with the discussion of sec. 4, which suggests that the error in the 
eigenvalue is proportional to ||Z|| ra , in which case n niin = cj log 10 ||Z||. We 
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have plotted this quantity as the dashed curve in fig. 4, choosing the constant 
c to be —5 and calculating || Z || with the Q-space eigenstate nearest to e 2 . 
As can be seen the dashed curve agrees very well with the exact result. 

6 Concluding Remarks 

We have derived an iteration method (GLS) for effective interactions as a 
generalization of the Lee-Suzuki method so that one can apply it to a system 
with arbitrary non-degenerate unperturbed energies. It has been proved 
that the iterative solution can be constructed in the framework of the Q- 
box formalism of Kuo et al. [3,5] in spite of the non- degeneracy of the 
unperturbed energies. The convergence in the GLS method depends on the 
overlaps of the exact P-space eigenstates and the unperturbed states. If the 
overlaps are small, the GLS approach does not always converge. 

The GLS method has further been generalized so that different projection 
operators can be used in each step of iteration. It has been shown that if we 
construct the projection operators using the P-space eigenstates determined 
self-consistently at each iteration step, the convergence is governed only by 
energy ratios of the differences between true and unperturbed energies and 
there is, in general, a unique way of distributing the true eigenvalues be- 
tween the sets E p and E q , where E p denotes the eigenvalues obtained from 
the P-space eigenvalue equation. This implies that, by appropriate choice 
of the unperturbed energies, the d eigenvalues obtained from the effective 
interaction can be tuned to generate any subset of the true eigenvalues. We 
have verified this theoretical prediction in the model calculations. We have 
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observed some divergent cases for the GLS, but the use of self-consistent 
projection operators in the SCGLS renders the iterations convergent. It also 
accelerates the convergence compared to the usual L-S method. 

This work was initiated during the authors' visit to the Nuclear Theory 
Institute at the University of Washington and we thank the Institute for 
the hospitality extended to us. Partial support from the US Department of 
Energy under contracts DE-FG02-88ER40388 and DE-FG02-87ER40328 is 
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Table 1. Convergence of the eigenvalues obtained with the SCGLS 
method for the model Hamiltonian with x = 0.05 









Number of iterations 




Exact 


ei 




2 


4 


6 


8 


10 


eigenvalues 


1.0 


2.0 


0.8915 


c 


c 


c 


c 


E 1 






2.2177 


2.2158 


c 


c 


c 


E2 


1.0 


3.0 


0.8966 


0.8904 


c 


c 


c 


E 1 






2.6998 


2.8553 


2.8619 


c 


c 


E 3 


1.0 


9.0 


0.9570 


0.8900 


0.8904 


c 


c 


E 1 






6.8904 


8.7779 


c 


c 


c 


£4 


2.0 


3.0 


0.8681 


0.5324 


2.5126 


2.2268 


2.2161* 


E2 






2.7209 


2.8541 


2.8466 


2.8621 


c 


E 3 


2.0 


9.0 


0.9677 


1.2027 


2.0943 


2.2106 


2.2155** 


E 2 






6.8689 


8.7827 


8.7816 


c 


c 


E 4 


3.0 


9.0 


1.6722 


3.6832 


2.8731 


2.8625 


c 


E3 






7.6619 


8.7861 


c 


c 


c 


E 4 



The notation c indicates convergence to four decimal places. The exact 
eigenvalues here are E 1 = 0.8905, E 2 = 2.2157, E 3 = 2.8622 and £ 4 = 8.7817 
and the overlaps of the exact eigenfunctions with the model space are 0.973, 
0.902, 0.121 and 0.004 respectively. 

Convergence to four decimal places is obtained after 12 iterations. 
** Convergence to four decimal places is obtained after 11 iterations. 
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Table 2. Convergence of the eigenvalues obtained with the SCGLS 
method for the model Hamiltonian with x = 0.20 







Number of iterations 




Exact 


ei 




2 


4 


6 


8 


eigenvalues 


0.0 


2.5 


-3.0221 


-0.2385 


c 


c 


E 1 






4.9933 


2.5150 


2.5793 


c 


E2 


0.0 


5.5 


-0.1608 


c 


c 


c 


E 1 






5.6553 


c 


c 


c 


E 3 


0.0 


9.0 


-0.1095 


-0.1495 


c 


c 


E 1 






8.8437 


c 


c 


c 


E 4 


2.5 


5.5 


2.0849 


2.5784 


c 


c 


E2 






5.4309 


c 


c 


c 


E 3 


2.5 


9.0 


2.5823 


c 


c 


c 


E 2 






8.9173 


c 


c 


c 


E 4 


5.5 


9.0 


1.4692 


6.3648 


5.6463 


c 


E 3 






8.9215 


8.9254 


c 


c 


E 4 



The notation c indicates convergence to four decimal places. The exact 
eigenvalues here are E l = -0.1496, E 2 = 2.5794, E 3 = 5.6451 and £ 4 = 
8.9251 and the overlap of the exact eigenf unctions with the model space are 
0.715, 0.295, 0.759 and 0.231 respectively. 
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Figure Captions 



Figure 1. The error in the calculated eigenvalues as a function of the number 
of iterations. Smooth curves are drawn through the points to guide the 
eye. The curves are labelled by the method employed to solve the model 
Hamiltonian problem with x = 0.05. 

Figure 2. As for fig. 1, but with x = 0.20. Here the calculations are 
converging to the first and third eigenvalues. 

Figure 3. As for fig. 1, but with x = 0.20. Here the calculations are 
converging to the second and third eigenvalues. 

Figure 4. Number of iterations required for an accuracy of one part in 10 4 
as a function of the unperturbed energy, e 2 , with x = 0.20. The full curve 
is obtained with the SCGLS method and the dashed curve gives the results 
expected from the convergence analysis, namely — 5/log 10 || Z ||. 
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